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Abstract 

The flow about slender, pointed bodies can be characterized by different 
states with angle of attack. At moderate-to-high angles of attack (qj ~ 40 ), 
a steady, asymmetric vortex pattern develops along the body, leading to a net 
lateral force. At higher angles of attack (a « 60°), the alt-end of the body devel- 
ops an unsteady von Karman shedding. As the angle of attack approaches 90 , 
the entire body length exhibits a time-dependent vortex shedding pattern. The 
current work uses three-dimensional, thin-layer Navier-Stokes simulations to 
investigate the physical mechanisms of asymmetric vortex shedding at a = 40° 
and a = 60°. The development of an asymmetric vortex pattern via a convec- 
tive instability mechanism is investigated using tip bumps, surface roughness, 
and tip curvature. It’s found that surface roughness simulations can incite an 
asymmetric vortex state at ci — 60° which is consistent with the application 
of a tip bump, and the experimentally observed flowfield. The unsteady von 
Karman vortex shedding on the aft portion of the body is also well resolved. 
The use of surface roughness did not incite a flow asymmetry at t* = 40 , and 
it was necessary to simulate tip curvature at this angle of attack in order to 
generate an asymmetric vortex state. 


1 Introduction 

It is well known that the wake behind a cylindrical body aligned perpendicular 
to the oncoming flow can develop an alternating pattern of strong vortices. This 
alternate shedding of vortices from the body results in an instantaneously asymmetric 
flowfield, while in the mean* the flowfield remains symmetric. The alternate shedding 
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is asserted in the United States under Title 17, U. S. Code. The U. S. Government has a royalty-free 
license to exercise all rights under the copyright claimed herein for Governmental purposes. All 
other rights are reserved by the copyright owner. 

’mean in this work refers to time-averaging. 
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of vortices from a circular cylinder is commonly referred to as von Karrnan vortex 
shedding, or a von Karrnan vortex street. If a slender, pointed body of circular cross- 
section, such as a cone-cylinder or ogive cylinder, is placed at an angle of attack 
relative to the oncoming flow, it’s possible to observe a steady asymmetric flow state 
relative to the lateral symmetry plane of the body. This introduces the counter- 
intuitive notion that a nominally axisvmmetric body can develop an asymmetric 
flow, not just instantaneously as in the flow about a circular cylinder in crossflow, 
but in a time-averaged sense*. 

The flow about a pointed, slender body of revolution at angle of attack can thus 
be broken down into several (experimentally-observed) [2-9] regimes (cf. Fig. T). The 
angles-of-attack where these flow regimes transition is dependent upon a number of 
factors, including body geometry, Mach number, Reynolds number, etc. At relatively 
low angles of attack, roughly a < 20°, the flowfield remains symmetric with respect 
to the lateral symmetry plane of the body, even if the flow is perturbed slightly. The 
flow in this regime may remain attached, or boundary layer crossflow separation can 
occur if the angle of attack is sufficient, leading to a pair of vortices on the leeward 
side of the body which are symmetric relative to the lateral symmetry plane. At the 
opposite angle of attack extreme, 70° < a < 90° (again approximately), the flow 
about the length of the body exhibits von Karrnan vortex shedding. In between these 
extremes, it’s possible for a nominally symmetric body of revolution to develop a non- 
zero mean lateral force. In the low-moderate angle-of-attack range, 20° < a < 40°, 
the variation of mean lateral force with roll angle forms a continuous periodic curve. 
In this regime, as the angle of attack is increased, the maximum absolute value of 
the lateral force also increases. As the angle of attack is increased into the moderate- 
high range, 40° < a < 70°, the variation of lateral force with roll angle can become 
nearly a square- wave, and the case of a symmetric mean flow becomes very difficult 
to achieve. This can be referred to as a “bi-stable” regime, or a bi-stable state; the 
lateral force becomes independent of roll angle except for a “switching 1 from positive 
to negative (and vice versa) as a critical roll angle is passed. 

Many of the previous experimental investigations noted that small changes near 
the tip of the model had a large effect on the observed vortex shedding pattern. 
Degani and Schiff [10], Degani [11], and Degani and Tobak [7-9] proposed that im- 
perfections near the tip of the model could incite a convective instability in the flow, 
leading to the development of steady asymmetric vortex patterns at high angles of 
attack. A convective instability is a sub-class of more general instability responses 
(cf. [12]). This is shown schematically in Fig. 2 as the response to an impulsive 
perturbation. The behavior in Fig. 2a is stable, in that the response decays in all 
directions with time. Fig. 2b shows an example of an absolutely unstable flow. Here 
the response to the initial disturbance grows in both the upstream (x/t < 0) and 
downstream ( x/t > 0) directions. A convectivelv unstable behavior shows growth 

*It is also possible for a nominally symmetric body to exhibit a steady asymmetric flow state in 
crossflow due to an asymmetric boundary layer transition to turbulence(cf. Ref. [1]). however the 
current work is only concerned with asymmetric vortex states. 




continuous with roii orientation 




unsteady, symmetric mean flow 
(von Karman shedding) 


Figure 1: Schematic of flow states about slender pointed bodies. The shaded lines 
represent vortices rotating in opposite directions. 

only in a downstream traveling packet. Given this response, if the perturbation is 
removed, the flow eventually reverts back to the state before the perturbation was 
applied. 

Computational Fluid Dynamics (CFD) simulations have been used to demonstrate 
the possibility of a convective instability mechanism for these class of flows. These 
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(a) Stable (b) Absolutely Unstable (c) Convectively Unstable 

Figure 2: Schematic of responses to an impulsive perturbation. 

investigations have primarily emphasized using the controlled environment of a simu- 
lation to perform numerical experiments that cannot be performed in a wind tunnel. 
By perturbing the geometry near the tip of the body, Degani et al. [11. 13] and Levy 
et al. [14] were able to demonstrate that the flowfield in the low-moderate angle-of- 
attack range behaves as if it possesses a convective instability. In the work of Degani 
et al., an asymmetric bump perturbation was applied to the tip of an ogive cylinder 
body until an asymmetric flow state developed, then removed, and the flow reverted 
to the symmetric state that existed before the perturbation was applied. Levy et 
al. noted that in the low-moderate angle of attack range increasing the size of the 
perturbation led to an increase in the flow asymmetry, consistent with the idea of 
a convective instability. Similar investigations in the moderate-high angle-of-attack 
regime [7] have indicated that the onset of von Karman vortex shedding in this regime 
may be due to an absolute flow instability, although the numerical results are more 
difficult to interpret in this angle-of-attack regime. Some researchers have been able 
to simulate an asymmetric mean flow in the moderate-high regime without perturbing 
the flowfield in any manner [15,16], however these simulations used conical flow, or 
other simplifying assumptions, which do not allow for the possibility of a convective 
instability. 

In Ref. 17, the author found that ogive cylinders placed at 40° angle of attack 
would exhibit an unsteady variation in normal force at Rep = 80.000 and 200,000, 
where D is the diameter of the cylindrical portion of the body. These variations were 
occurring at harmonics of the Strouhal frequency predicted by the Independence 
Principle, i.e. St n = Sf 2 o sin or, and hence were incorrectly related to a possible von 
Karman shedding mechanism. Subsequent investigations [18] of crossflow separation 
in two and three-dimensions, and at several Reynolds numbers using varying grid 
densities, showed that these unsteady variations were caused by an instability in 
the free shear layer shed from the body. In a physical flow, these free shear layer 
instabilities would eventually transition and lead to turbulence in the wake behind 
the body, as examined experimentally by Roshko [19] and later bv Bloor [20]. In 
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Ref. 17, it was also noted that the amplification of round-off errors was sufficient to 
lead to a von Kantian vortex shedding along the length an ogive cylinder at a = 60° 
and Rr D = 200, 000, and no geometric or other How perturbations was required to 
initiate this state. 

The current work uses C’FD simulations to investigate the physics of vortex shed- 
ding about pointed slender bodies; both the steady, asymmetric vortex patterns at 
moderate angles of attack, and the unsteady, von Karman vortex shedding that oc- 
curs at high angles of attack. Since a large database of experimental results exists for 
steady asymmetric vortex shedding for Reynolds number Reo = 10, 000 — 200, 000, 
the current work is performed in this same Reynolds number range. While it has been 
established that imperfections near the tip of a pointed body can incite a convective 
instability in the flowfield leading to an asymmetric vortex pattern, the physical mech- 
anisms by which this convective instability affects the changes on the flowfield are still 
not well understood. Without a detailed understanding of the physical mechanisms 
involved, designers seeking to eliminate or harness the lateral control power of asym- 
metric vortex shedding are left to test and develop concepts based solely on physical 
intuition. The current work seeks to provide insight into the physical mechanisms in- 
volved by performing controlled numerical experiments comparing different methods 
of perturbing the flow; tip bumps, surface roughness, and tip curvature. 

2 Approach 

2.1 Numerical Algorithm 

This work deals with numerical modeling of two- and three-dimensional separated 
flowfields. In order to accurately simulate the physics of these flows, it is necessary 
to use the Navier-Stokes equations as a physical model. The complete Navier-Stokes 
equations are often simplified using the “thin-layer” approximation [21] when perform- 
ing high Reynolds number flow simulations. Neglecting the crossflow contributions 
to the viscous stress tensor still provides satisfactory results for the computed mean 
drag, and von Karman shedding frequency at the flow conditions of interest in this 
work (cf. Ref. 18), and all of the simulations in the current work utilize this simplifi- 
cation to reduce the computational costs. At the Reynolds number considered in the 
current work, the boundary layer remains laminar, and a numerical turbulence model 
is not required. 

In order to accurately resolve the viscous stress terms, a fine grid spacing is re- 
quired in each coordinate direction that the viscous stresses cannot be neglected. 
This grid spacing makes the use of implicit methods desirable when using finite- 
difference algorithms because of the stiffness of the Navier-Stokes equations. The 
implicit schemes used in this work are all based on performing a time-linearization 
of the non-linear flux vectors, and obtaining a so-called “delta formulation” of the 
governing equations [22]. The backward Euler implicit time-integration scheme is 



used throughout.. This is a first-order accurate scheme in time, however the allowable 
numerical time-step due to stability considerations for these flows is much smaller 
than required to resolve the physics of the unsteady von Karman shedding, and a 
first-order scheme has been found to provide satisfactory results (cf. Refs. 17. 18). 

In order to numerically integrate the governing equations, the spatial derivatives 
must be replaced by suitable finite-difference formulas. A straight-forward approach 
is to replace the inviscid derivatives with second-order, central-difference operators, 
and the viscous derivatives with centered, second-order mid-point operators. Unfortu- 
nately, this produces a scheme which is moderately unstable in three dimensions [23]. 
This weak instability can be overcome however with the addition of numerical dissipa- 
tion of an order that doesn’t interfere with the accuracy of any physical mechanisms. 
The explicit dissipation is scaled by the local flow properties to reduce it’s influence 
within the boundary layer, where the natural viscous dissipation of the Navier-Stokes 
equations is suitable. 

This combination of central-differenced spatial terms and higher-order numerical 
dissipation leads to what is referred to here as a block form of the Beam and Warming 
algorithm [22,23]. The block formulation is preferred over the diagonalized algorithm 
of Pulliam and Chausee [24] due to an inherent asymmetry of the diagonalization pro- 
cess that can lead to spurious flow asymmetries for high-angle-of-attack simulations 
(cf. Ref. 25). 

2.2 Computational Grids 

Computations were performed about a tangent-ogive cylinder configuration with 
a 3.5 diameter forebody attached to a 7.0 diameter cylindrical aft section. This is the 
same configuration examined computationally by Degani et al. [7-9, 11], in previous 
work by the author [17], and experimentally by Lamont [5,6]. In previous simulations 
of vortex shedding about a circular cylinder in crossflow and a tangent-ogive cylinder 
at a = 60° and a = 80°, it was found that the von Karman vortex shedding frequency 
was not well predicted at Reo = 80,000 or Reo = 200,000 (cf. Ref. 17). Grid 
refinement studies [18] determined a minimum crossflow resolution (A0 = 1.0°, where 
(j) is the circumferential spacing) required in the separated flow region in order to 
predict the von Karman shedding frequency, mean drag, and crossflow separation 
location at these flow conditions. This recommended crossflow resolution is utilized 
in this work. This leads to a grid with a crossflow plane having 249 circumferential 
and 58 radial points for a full-body (periodic) configuration, and 69 points in the 
axial direction. A typical axial grid plane is shown in Fig. 3 (every third grid line is 
shown for clarity). The outer boundary of the computational domain is located fifteen 
body diameters from the surface. At the outer boundary a characteristic boundary 
condition is applied, while at the downstream exit plane a zeroth-order extrapolation 
of the flow quantities is performed, which is consistent with the highly convective 
nature of the flowfield. 
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Figure 3: Lateral symmetry plane of the tangent-ogive cylinder computational grid. 

3 Numerical Results 
3.1 Tip Bump 

As was described in Sec. 1, the flow about a pointed slender body at moderate-to- 
high angles of attack can develop a steady asymmetric vortex pattern which produces 
a non-zero mean lateral force on the body. Previous CFD studies [11,13,14], were 
able to simulate an asymmetric vortex pattern on a slender body by placing a bump 
on the body very close to the tip of the nose. In [17], a similar numerical simulation 
was performed at a = 60° and R&o = 200,000. This simulation utilized a similar 
crossflow resolution (A 0 = 3.0°) as the previous studies by Degani et al. and Levy 
et al. The tip-bump simulation from [17] will be used as a baseline for the current 
work, and hence is first briefly described. 

The computed time-histories of the lateral and normal forces are plotted in Fig. 4. 
Initially, the geometry was unperturbed, and a von Karman vortex shedding pattern 
develops along the length of the body. This can clearly be seen in the variation 
of lateral force. The computed vortex shedding frequency and the circumferential 
location of the crossflow separation are not in good agreement with experimental 
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investigations of shedding from a two-dimensional cylinder (cl. Ret. 18). At t 266, 
where i is the non-dimensional time, a hump of height h max = 0.0 Id. where d is the 
local diameter, was placed at r/D = 0.06 and 90° from the windward symmetry plane 
(cf. Appendix A). After the simulation had progressed through several cycles, the 
perturbation was removed at l = 424. 



(a) Lateral force (b) Normal force 

Figure 4: Tip bump perturbation force coefficient histories = 0.20, a = 60°, 

Re D = 200,000). 

When the bump is placed on the nose the flow develops a non-zero mean lateral 
force, while still maintaining vortex shedding about the aft-end of the body. The 
presence of the bump also causes the mean normal force to increase. This is consistent 
with experimental observations (cf. Ref. 4), which note that the greater the lateral 
asymmetry, the greater the normal force. The computed flowfield without the bump 
present contains two main vortices on the leeward side of the body that convect at a 
small angle to the longitudinal axis of the body and alternate in time. Figure 5 shows 
an instantaneous snapshot of the unsteady streaklines [26] after the tip perturbation 
has been applied. In the presence of the bump, the tip vortex on the perturbation 
side (R) quickly curves away from the body, while the vortex on the opposite side (L) 
remains attached for a short distance, then curves away from the body. This steady, 
asymmetric pattern of vortices produces the non-zero mean lateral force. This steady 
vortex pattern is consistent with those observed experimentally [7,27] and previously 
simulated numerically by Degani [13] using a similar bump perturbation. Downstream 
of these two steady, asymmetric vortices the flowfield remains unsteady with the tip 
bump in place. 

When the bump is removed, the flowfield reverts to a symmetric time-averaged 
flow. This is consistent with the presence of a convective instability mechanism, as 
has been proposed by Degani and Tobak [7-9]. 
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Figure 5: Instantaneous snapshot (t = 423) of streaklines with tip bump present 
(A/*, = 0.20, a = 60°, Re D = 200, 000). 

3.2 Surface Roughness 

The placement of a fixed bump on the tip of the ogive cylinder is a useful tool 
for examining the response of the flowfield to a particular perturbation, however it 
does have limitations. A criticism of the method is that it places a relatively large 
perturbation (compared to the local radius of curvature) where the flowfield is known 
to be most sensitive. At low angles of attack it is necessary to use large bump heights 
to produce the same level of lateral force as observed experimentally (cf. [14]). 

The use of a large surface perturbation to “trigger” the convective instability 
mechanism may effect this change via a different physical mechanism than occurs 
experimentally. The large perturbation can cause the boundary layer to separate at 
the bump, and become entrained in the outer, inviscid flowfield. This could bias the 
flow towards a crossflow separation driven mechanism for asymmetric vortex shedding. 

In the current work, it was investigated whether surface roughness can produce the 
same asymmetric flow states as observed in both previous numerical simulations with 
fixed tip bumps, and experimentally. The use of surface roughness more accurately 
mimics the physical geometry encountered with an experimental model, and is easier 
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to simulate numerically. The modeling of a bump near the tip can cause' problems tor 
both grid generation and numerical stability, while small surface perturbations are 
easier to generate and do not appreciably skew the existing grid cells. 

The flow about a tangent-ogive cylinder at <v = 60° was computed using a crossflow 
grid resolution of A</> — 1.0°. The lateral and normal force histories for the simulation 
are plotted in Fig. 6. The simulation initially had a “clean’" geometry, and again the 
flowfield develops a von Karman vortex shedding before the perturbation is applied. 
The variation of lateral force (and normal force) is larger than in the previous tip 
bump simulation, which is indicative of stronger vortices, and the vortex shedding 
frequency is lower. Both of these observations are consistent with the increase in 
crossflow grid resolution in the current surface roughness simulation compared to 
the prior tip bump studies. The vortex shedding frequency follows very closely the 
curve St a = St 2 D sin a suggested by the Independence Principle. At 1 = 155.5 
a random distribution of surface roughness with A max = 1.0 x 10 4 was applied 
(cf. Appendix A). After the surface roughness is applied, the flowfield develops a 
non-zero mean lateral force in the positive, or pilots-right, direction. The sign of 
the lateral force (configuration of the asymmetric vortex pattern) would essentially 
be chosen at random with this method. At t = 228.5 the surface roughness was 
removed, and the flowfield reverted back to a mean symmetric flow, again consistent 
with the behavior of a convective instability. 



(a) Lateral force (b) Normal force 

Figure 6: Surface roughness force coefficient histories (Mx> = 0.20, a = 60°, Reo = 
200 , 000 ). 

Figure 7 shows an instantaneous snapshot of the unsteady streaklines computed 
with surface roughness applied. As in the bump simulation, one tip vortex (L) lifts 
rapidly from the body, while the lifting of the tip vortex on the opposite side (R) is 
delayed. An interface between the space-fixed oblique shedding and the unsteady von 
Karman vortex shedding on the aft-end of the body can be seen. 

Another useful tool for examining the asymmetric flowfield, and it’s response to 
perturbations, is the variation of sectional lateral force with both axial distance along 
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Figure 7: Instantaneous snapshot (t — 205) of streaklines with surface roughness 
present (Moo — 0.20, a = 60°, Rep = 200,000). 

the body and non-dimensional time. This data is presented as a three-dimensional 
carpet plot with sectional lateral force being the dependent variable. The sectional 
lateral force variation is plotted for both the tip bump simulation and the surface 
roughness simulation in Figs. 8 and 9 respectively*. The lateral force maps are similar 
for the two simulations. In both cases the unsteady vortex shedding is clearly visible 
on the aft-end of the body (x/D > 6). The frequency of the vortex shedding is 
higher, and the amplitude of the variation lower, in the tip bump simulation, which 
is consistent with its lower crossflow grid resolution. On the tangent portion of the 
nose, both simulations show a steady lateral force which originates as a perturbation 
from the nose tip. In the tip bump simulation the lateral force is negative (towards 
the pilots left), while in the surface roughness simulation the sign is reversed. Both 
simulations also show a high amplitude, high frequency variation of lateral force near 
the nose-cylinder junction ( x/D — 3.5). In the surface roughness simulation this 
variation is difficult to see because it’s superimposed on the steady lateral force on 
the nose region. These high frequency peaks may be related to the onset of vortex 
shedding on the cylindrical portion, and the high frequency shear-layer instabilities 
discussed in Sec. 1. 

‘Note that Figs. 8 and 9 do not have the same time scale. Refer to Figs. 4 and 6 to determine 
the perturbation application for the respective simulations. 
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It, isn't possible to see in the resolution of Figs. 8 and 9. but very close to the 
tip point., the fiowfield is highly unsteady in the tip bump simulation, while in the 
surface roughness simulation the flow is steady (and asymmetric). This unsteadiness 
is caused by the tip bump forcing a crossflow separation. Figure 10 contains the 
helicity density contours in a crossflow plane located at x/D — 0.06 for both the tip 
bump and surface roughness simulations. As can be seen, the tip bump causes an 
irregular, highly unsteady configuration of the main vortices, while the asymmetric 
vortices in the surface roughness simulation remain steady. The irregularities in the 
free shear layer due to the shear-layer instabilities in the surface roughness simulation 
are also visible. The steady asymmetric flow pattern that is observed on the entire 
forebody in the surface roughness simulation is an indication that the flow has been 
perturbed (slightly) into a “stable” flow state, not one that is artificially induced. 

When the steady asymmetric vortex pattern and the unsteady von Karman shed- 
ding exist simultaneously, it’s hypothosized that a front separates the two states, 
possibly similar to the weak pressure fronts described by Williamson [28]. As seen in 
the numerical simulations, and through physical arguments (cf. [8]), this front must 
remain stationary relative to the body (Fig. 11). If the tip perturbation is removed, 
the new flow state can also be interpreted in the light of a response to a perturbation, 
in this case removing the tip imperfections. Since the von Karman shedding appears 
as an absolute instability, the response moves faster than the body travels, hence the 
pressure front should move upstream until the entire body experiences unsteady von 
Karman shedding. 

3.3 Tip Displacement 

It was attempted to generate an asymmetric flow at a = 40° without resorting 
to the use of a tip bump. It was first attempted to utilize the surface roughness 
simulation that was effective at a = 60°, as described above. Initially, the same 
roughness height of A max — 1.0 x 10 4 was applied, however this did not generate 
an appreciable lateral force. The roughness height was increased to A max = 0.001 
everywhere, yet still it was not possible to generate an asymmetric flow state. Finally, 
a configuration was created which maintained a roughness height of A max = 0.01 on 
the nose of the body (x/D < 1.0) and A mttX = 0.001 elsewhere, yet this still did not 
generate an appreciable lateral force on the body at this angle of attack. 

Another geometry modification which is consistent with experimental models is an 
asymmetry of the sharp tip. Moskovitz et al. [29] reported that the tip of their model 
appeared curved when viewed under a microscope. A configuration was generated 
where the tip point was offset 0.001 D in the y-direction from the lateral symmetry 
plane, while the remainder of the geometry was unmodified. The flow about this 
geometry was simulated at ci = 40°. The lateral force immediately began to diverge 
from zero after the tip was moved, and eventually the flow reached a steady asymmet- 
ric state. A similar effect was found at a - 60°, however a comparable computation 
to the surface roughness simulation described above has not been completed to date. 
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Figure 10: Helicitv-density contours at x/D = 0.06 (A/oo — 0.20, a — 60°, Reo — 

200 , 000 ). 


The effect of the tip displacement on the the lateral force evolution can be traced to 
the perturbation response being convected down the length of the body. Figure 12 
shows the variation of sectional lateral force along the body with time. The response 
convects downstream creating an asymmetric vortex strength, with the asymmetry in 
the positive direction for 3 < x/D < 6 and in the negative direction for 6 < x/D < 10. 
The final distribution of sectional lateral force is in relatively good agreement with 
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Figure 11: Schematic of front separating steady and unsteady flow regions. 


the experimentally observed modes (cf. Refs. 4,6,14). 

From the off-surface streamline pattern (cf. Fig. 13) it was noted that the tip dis- 
placement in this configuration did not cause an asymmetric oblique vortex shedding 
pattern to develop, rather the two tip vortices have different strengths on opposite 
sides of the body. This is a prelude to the asymmetric vortex pattern that can de- 
velop, and is consistent with observations of a continuous “sine-wave distribution” 
of lateral force with roll angle. It was attempted to cause the vortices to “lift-off” 
the body by increasing the tip displacement. This did lead to a larger flow asymme- 
try, but the vortices remained nearly parallel to the body longitudinal axis. Similar 
results were obtained by combining a tip curvature simulation with a surface rough- 
ness of A max = 0.001 everywhere, and combining tip curvature with a slight sideslip 
angle (0 = 2.0°). The trends observed (that a larger perturbation leads to a larger 
asymmetry) are consistent with previous numerical results from Levy et al. [14] using 
tip bumps to perturb the geometry. Levy et al. were able to cause the vortices to 
lift off the body at this angle of attack, however in that case the flowfield was fully 
turbulent, which may affect the behavior. Similar to the results obtained using sur- 
face roughness simulations, utilizing a tip displacement did not cause the crossflow 
separation location near the tip of the body to become unsteady. 


4 Summary 

Simulations of asymmetric vortex states about a tangent-ogive cylinder configu- 
ration at high angles of attack were performed using three means of geometrically 
perturbing the body; tip bumps, surface roughness, and tip curvature. It is seen 
that the tip bump does provide a useful tool for simulating the cumulative effects 
of tip asymmetries, however it can artificially excite the asymmetric state, especially 
near the tip of the nose. A method of simulating surface roughness was found to 
produce a steady asymmetric state on the forebody of the ogive cylinder at a = 60°, 
with unsteady von Karman shedding on the aft end of the body. At a = 40°, the 
surface roughness was ineffective in producing an asymmetric vortex configuration, 
however a method of simulating tip curvature did produce a pair of leeward vortices 
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Figure 12: Variation of lateral force with tip displacement present (Moo = 0.20, 
a = 40°, Re D = 200,000). 



Figure 13: Off-surface streamlines with tip displacement present (Moo = 0.20, a - 
40°, Re D = 200,000). 

with different strengths. Both of these methods are considered non-intrusive , in 
that any changes to flowfield near the tip are induced rather than forced, and the 
flowfield remains steady. However, it was not possible to simulate large flow asym- 
metries at a = 40° with the tip curvature method. The use of the surface roughness 
or tip curvature methods are preferred over the tip bumps as they are simpler to 
implement numerically, and can be used for detailed investigations into the physical 
mechanisms causing the vortex asymmetries without forcing the crossflow separation 
to be displaced. 
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A Geometry Perturbations 

A.l Bump Geometry Definition 




(a) Front view (y-z plane) (b) Top view ( x-y plane) 

Figure 14: Bump geometry. 

The bump geometry is defined as a circular perturbation following Fig. 14. h max 
is the maximum bump height, and is placed at an axial location x b . The x-axis is 
assumed to coincide with the longitudinal axis of the body. The magnitude of the 
bump radius, |r&|, is also specified. A value of |r ft | = 0.5|r(x)| has been found to give 
good results for tip bumps. The longitudinal distribution of bump contours is created 
by defining a bump height distribution in x, and building planes of circular bumps, 



Where l b is the half-length of the bump perturbation. A value of £ b = 2|r(x fe )| has 
given good results with tip bumps. This leaves the origin of the bump circle to be 
solved from 

\r(x) \ + h(x) = |r„| + \r b \ 

In the current work the bump was always placed 90° from the windward symmetry 
plane. 

In this manner the existing grid points (x) in the crossflow plane will be modified 
only if they are located within the region of the bump that protrudes from the body, 
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i.e. 


if (|x — r „ | < |r*&| ) then 
x = r„ 4 - r h 
fi 


A. 2 Surface Roughness Simulation 

The effects of surface roughness were modeled by randomly applying small per- 
turbations to the geometry surface, 


r' = r(x)(l + A) — A max < A < A max 

The maximum possible deviation between successive points is thus 2A max . The per- 
turbation was only applied to half the surface points, 


A = 


0 

±A max rand[0, 1] 


if rand[0, 1] < 0.5 
otherwise 


where rand [0,1] is a pseudo-random number between 0 and 1 inclusive. The pertur- 
bation is also chosen to be positive or negative at random. The random numbers in 
the above equation (including the ± test) are considered to be distinct. 
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